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REMOTE RESERVOIR RESISTIVITY 
MAPPING 

This application claims the benefit of U.S. provisional 
application No. 60/154,114 filed on Sep. 15, 1999. 

FIELD OF THE INVENTION 

This invention relates generally to the field of geophysical 
prospecting. More particularly, the invention relates to sur- 
face measurement of subsurface geologic formation electri- 
cal resistivity. Specifically, the invention is a method of 
combining seismic and electromagnetic data to prospect for 
subsurface formations that contain hydrocarbons. 

BACKGROUND OF THE INVENTION 

Remote mapping and analysis from the surface of the 
earth of hydrocarbons reservoired at depth remains a diffi- 
cult technical task. This is so despite recent advances in 3D 
seismic imaging, seismic direct hydrocarbon indicator 
(DHI) and amplitude variation with offset (AVO) analyses, 
and seismic attribute mapping and interpretation. Seismic 
detection difficulties arise in part from the fact that the 
mechanical properties of reservoirs, to which the seismic 
probe responds, are often only slightly modified when 
hydrocarbons replace formation waters, especially in the 
case of oil. The modification may be of the order of only 
10's of percent. Subtle mechanical effects related to seismic 
wave propagation and reflection can mask DHI and AVO 
signatures or even produce misleading signatures. For 
example, low gas saturation in water sands can produce false 
seismic DHIs. Because of such effects, drill-well success 
rates are too low and exploration costs are too high in many 
basins. In addition, rapid and low-cost assessment of dis- 
covered undeveloped hydrocarbon reserves requires good 
knowledge of reservoir properties at large distances from the 
discovery well. Acquiring this knowledge is problematic 
using only seismic data. There is an urgent need to remotely 
measure and map other reservoir formation properties that 
are sensitive to hydrocarbons, and to combine interpretation 
of these other properties with interpretations of seismic data 
and their mapped attributes. One particularly important 
formation property is electrical resistivity, which is strongly 
related to the pore fluid type and saturation. 

The bulk electrical resistivity of reservoirs is often 
increased substantially when hydrocarbons are present. The 
increase can be of the order of 100's to 1000's of percent. 
However, increased formation resistivity alone may not 
uniquely indicate hydrocarbons. For instance, carbonates, 
volcanics, and coals can also be highly resistive. 
Nevertheless, spatial correlation of high formation resistiv- 
ity with potential traps imaged by seismic data, or with 
seismic DHI or AVO effects at reservoir depth, provides 
strong evidence of the presence of oil or gas and valuable 
information on their concentrations. For example, a low gas 
saturation high-porosity sandstone reservoir encased in 
shale can produce a strong seismic DHI and an AVO curve 
indicative of gas. However, it would also have low electrical 
resistivity and hence would be a high-risk drill-well pros- 
pect. 

Most hydrocarbon reservoirs are inter-bedded with shale 
stringers or other non-permeable intervals and hence are 
electrically anisotropic at the macroscopic scale. Thus, it is 
important to measure both the vertical (transverse) and 
horizontal (longitudinal) electrical resistivities of the reser- 
voir interval. Remote measurement of the vertical and 
horizontal resistivities of the reservoir interval, combined 
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with estimation of the resistivity of the non-permeable 
bedding, would provide quantitative bounds on the reser- 
voir's fluid content, such as the hydrocarbon pore volume. 
However, there is no existing technology for remotely 

5 measuring reservoir formation resistivity from the land 
surface or the seafloor at the vertical resolution required in 
hydrocarbon exploration and production. Based on the 
thicknesses of known reservoirs and predicted future needs, 
this required resolution would be equal to or less than two 

10 percent of depth from the earth's surface or seafloor. For 
example, this would resolve a 200-ft net reservoir thickness 
(vertical sum of hydrocarbon bearing rock thicknesses 
within the reservoir interval) or less at a typical 10,000-ft 
reservoir depth. 

15 Overviews of electromagnetic imaging technology are 
given by M. N. Nabighian (ed.), Electromagnetic Methods in 
Applied Geophysics, Vols. 1 & 2, SEG Investigations in 
Geophysics No. 3, 1988; A. G. Nekut and B. R. Spies, 
Proceedings IEEE, v. 77, 338-362, 1989; and by M. S. 

20 Zhdanov and G. V. Keller, The Geoelectrical Methods in 
Geophysical Exploration, Elsevier, 1994. Imaging of elec- 
trically conductive objects such as ore bodies has been the 
dominant application for electromagnetic methods. In appli- 
cations for hydrocarbon exploration, most of the technology 

25 was developed to image large geological structures in 
regions where seismic data are low in quality or are absent, 
and little other geological or geophysical information is 
available. 

Direct exploration for hydrocarbons using surface-based 

30 electromagnetic imaging has been attempted since the 
1930s, but with little commercial success. This lack of 
success is due to the low spatial resolution and the ambigu- 
ous interpretation results of current electromagnetic 
methods, when applied in stand-alone and spatially under- 

35 sampled ways to the geological imaging problem. Low 
subsurface resolution is one consequence of the diffusive 
nature of the low frequency electromagnetic waves, that is, 
below 1 kHz, required to penetrate the earth to reservoir 
depths. The vertical resolution of such electromagnetic 

40 waves is relatively insensitive to bandwidth, unlike the 
seismic case, but is very sensitive to the accuracy and 
precision of phase and amplitude measurements and to the 
inclusion of constraints from other data. That is, the uncon- 
strained geophysical electromagnetic data inverse problem 

45 is mathematically ill posed, with many possible geologic 
structures fitting electromagnetic data equally well. 
Consequently, the vertical resolution of unconstrained elec- 
tromagnetic imaging is typically no better than 10 percent of 
depth. This gives a resolution of only a 1000-ft net reservoir 

50 thickness at a typical 10,000-ft reservoir depth. However, 
within a given resolved layer, conventional resistivity mea- 
surement accuracy can be within a factor of two, which is 
adequate for oil and gas exploration. 
Electromagnetic technology that is applicable to direct 

55 reservoir imaging uses electrically grounded controlled 
sources to produce vertical and horizontal current flow in the 
subsurface at the reservoir depth. The five embodiments of 
this technology, well known within the electromagnetic 
imaging community, are: (1) the LOTEM method described 

60 by K. M. Strack, Exploration with Deep Transient 
Electromagnetics, Elsevier, 1992; (2) the SIROTEM 
method, described by Buselli in U.S. Pat. No 4,247,821; (3) 
CGG's TRANSIEL® system, described in U.S. Pat. No. 
4,535,5293; (4) the EMI method, described by Tasci et al. in 

65 U.S. Pat. No. 5,563,513; and (5) the WEGA-D method 
described by B. W. Smith and J. Dzwinel in WEGA-D 
SYSTEM®, WEGA-D Geophysical Research Ltd., 1984. A 
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newer version of WEGA-D named PowerProbe® has been arrays in their calculations. Maurer and Boerner discuss the 

developed by the Canadian company Enertec, a successor to more general problem of optimization of surface electro- 

WEGA-D Geophysical Research. All five methods suffer magnetic surveys for imaging subsurface targets, but do not 

from the vertical resolution limitation of approximately 10% discuss unipole, multiple radial bipole, or concentric ring 

of depth cited above, which makes them unsuitable for direct 5 dipole sources. 

reservoir imaging except for unusually thick reservoirs. This Conventional geophysical electromagnetic data process- 
resolution limitation results from one or more of the fol- ing finds the minimum earth structure, that is, the simplest 
lowing deficiencies in each method: (1) lack of means to resistivity model, which is consistent with the measured data 
focus the electromagnetic input energy at the target reser- within the experimental error bounds, but without explicit 
voir; (2) spatial undersampling of the surface electromag- 10 incorporation of a priori information. Incorporation of hard 
netic response fields; (3) measurement of only a few com- constraints into the data processing significantly improves 
ponents (usually one) of the multi-component spatial resolution and resistivity accuracy, which are not 
electromagnetic surface fields that comprise full tensor simply related to signal wavelength or bandwidth as in the 
electromagnetic responses at each receiver (except for seismic case. Examination of well log and other data shows 
WEGA-D/PowerProbe); (4) data processing using 1-D, 2-D, 15 that, in most cases, major seismic boundaries are also major 
or pattern recognition algorithms rather than full 3-D imag- resistivity boundaries. In addition, interpretation of seismic, 
ing methods; and (5) lack or paucity of explicit depth gravity, and magnetic data would provide good knowledge 
information and resistivity parameter values incorporated of the major lithologies present in a prospective area before 
into the data processing to constrain the inversion results. drilling. Applying constraints for a large number, (10's to 
Another serious limitation in these five methods is their 2 o 100's) of layers and other major geologic boundaries (for 
use of high-impedance contact electrodes and connecting instance, faults) would be novel for electromagnetic imaging 
wires, with greater than 1 Ohm total series resistance, to of hydrocarbon reservoirs. 

transmit the source current into the subsurface. This output Two previous methods have described the incorporation 
impedance is primarily a result of the small surface area of of seismic constraints to improve spatial resolution in low- 
the electrodes that contact (i.e. ground to) the earth. High 2 s frequency electromagnetic geophysical inversion. Although 
output impedance severely limits the electrical current at the not applied to hydrocarbon reservoir imaging, a method was 
reservoir depth, which in turn reduces the strengths of the developed by G. M. Hoversten et. al., Geophysics, v. 63, 
surface electromagnetic responses to the subsurface reser- 826-840, 1998a; and SEG Annual Meeting Expanded 
voir for a given source power. Current limitation due to Abstracts, v. 1, 425-428, 1998b to improve 2-D natural- 
high-impedance sources also results in reduced depths of 30 source electromagnetic (magnetotelluric) imaging of the 
exploration, especially in electrically conductive sedimen- base of salt structures in the offshore Gulf of Mexico, 
tary basins. The effective depth of electromagnetic explora- Vertical resolution of the salt base improves by a factor of 2 
tion increases as a fractional power of source strength, to 3 when the depth to the top of salt is constrained by 3-D 
between M 1/5 and M 1/3 for grounded electric dipole sources seismic data and when the salt resistivity is fixed. Natural- 
where M is the dipole moment, that is, current multiplied by 35 source methods such as that of Hoversten et al. lack the 
dipole length. The exponent depends upon which surface vertical resolution required for direct imaging of resistive 
field component is measured, but in general short -offset (or hydrocarbon reservoirs, because they measure the earth's 
"near-field") electromagnetic receiver responses have the response to the flow of horizontal subsurface electrical 
best sensitivity to deep targets, as shown in B. R. Spies, currents that are insensitive to regions of increased resistiv- 
Geophysics v. 54, 872-888, 1989. 40 ity. D. L. Alumbaugh and G. A. Newman, Geophys. J. Int., 
V. S. Mogilatov and B. Balashov,7. Appl Geophys., v. 36, v. 128, 355-363, 1997; and SEG Annual Meeting Expanded 
31-41, 1996; and Mogilatov's Russian patent 2,084,929-Cl Abstracts, v. 1, 448-451, 1998 have described the use of 
describe the use of surface electric concentric ring dipoles seismic constraints to improve resolution in cross-well elec- 
and radial electric bipoles. A. G. Tarkov, Bull Acad. ScL tromagnetic imaging within hydrocarbon reservoirs, in a 
U.S.S.R., Geophys. Ser., no. 8, 11, 1957, R. N. Gupta and P. 45 manner similar to that of Hoversten et al. for surface 
K. Bhattacharya, Geophysics, v. 28, 608-616, 1963, and by magnetotelluric data. However, the cross-well method 
A. Dey et al., Geophysics, v. 40, 630-640, 1975 describe the requires the existence of at least two wells that penetrate the 
use of opposite-polarity collinear surface electric bipoles reservoir. 

("unipoles"). However, ring electrodes described by Mogj- Estimation of the reservoir's fluid type, saturation, and 

latov and Balashov do not contain discussions of, much less 50 shaliness factor from surface geophysical measurements has 

calculations for, the optimum electrode dimensions needed been previously conducted using only seismic reflection 

to maximize the vertical electric field or current density at data, in particular various seismic interval attributes 

the target (reservoir) depth. The unipole methods described (amplitude widths, ratios, phases, etc.). Here, the shaliness 

by Tarkov, Gupta, Bhattacharya, and Dey et al. do not factor is the ratio of net hydrocarbon bearing zone thickness 

include discussions of or calculations for the effects of 55 (pay) to gross reservoir thickness. It is well known in the 

changing the source frequency, or the effects of using industry that the electromagnetic response of a vertically 

finite-length unipoles (second electrodes not at infinite layered earth depends on the direction of the resistivity 

distance), on the optimum configuration needed to maximize measurement. See, for instance, M. S. Zhdanov and G. V. 

the vertical electric field or current density at the target KeUer (1994, op. cit.). However, there is no existing remote 

depth. 60 (surf ace -based) electromagnetic method for measuring both 

S. K. Verma and S. P. Sharma, Geophysics, v. 60, the separate vertical and horizontal resistivities of a reservoir 

381-389, 1995 and H. Maurer and D. E. Boerner, Geophys. interval at depth. Directional resistivity measurements for 

J.Int.,v. 132,458-468, 1998 discuss optimization of surface reservoirs have been restricted to in-situ methods, such as 

electromagnetic source array configurations in order to best well logging. 

focus energy onto subsurface targets. However, Verma and 65 Specific technologies for indirect electromagnetic detec- 

Sharma restrict their discussion to subsurface conducting tion of reservoired hydrocarbons at depth have also been 

layers, and do not include unipole or concentric ring dipole developed, but these rely on the detection of electrically 
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altered zones ("chimneys") above reservoirs caused by the 
purported slow leakage of hydrocarbons upward from the 
reservoir. The existence and relationships of alteration chim- 
neys to reservoired hydrocarbons have not been unequivo- 
cally demonstrated. Changes in resistivity (increases and 
decreases) and polarizability (or induced polarization) are 
claimed by the practitioners of chimney detection to occur at 
various locations within such chimneys. Electromagnetic 
methods to locate chimneys were developed by Sternberg et 
al., as described in their U.S. Pat. No. 4,446,434, and Tasci 
et al, as described in their U.S. Pat. No. 5,563,513. The 
TRANSEEL® and WEGA-D/PowerProbe systems can also 
be used to detect hydrocarbon chimneys. These methods 
suffer the same depth resolution limitations as listed above, 
for the reasons cited in the preceding paragraph. 

SUMMARY OF THE INVENTION 

The invention is a method for surface estimation of 
reservoir properties of a subsurface geologic formation. 
First, the location of and average earth resistivities above, 
below, and horizontally adjacent to the subsurface geologic 
formation are determined, using geological and geophysical 
data in the vicinity of the subsurface geologic formation. 
Second, the dimensions and probing frequency for an elec- 
tromagnetic source are determined to substantially maxi- 
mize transmitted vertical and horizontal electric currents at 
the subsurface geologic formation, using the location and 
the average earth resistivities. Next, the electromagnetic 
source is activated at or near surface, approximately cen- 
tered above the subsurface geologic formation and a plural- 
ity of components of electromagnetic response are measured 
with a receiver array. Next, geometrical and electrical 
parameter constraints are determined, using the geological 
and geophysical data. Finally, the electromagnetic response 
is processed using the geometrical and electrical parameter 
constraints to produce inverted vertical and horizontal resis- 
tivity depth images. 

In an alternative embodiment, the inverted resistivity 
depth images may be further combined with the geological 
and geophysical data to estimate the reservoir fluid and 
shaliness properties. 

In a further alternative embodiment, the average earth 
resistivities above, below, and horizontally adjacent to the 
subsurface geologic formation are verified using the plural- 
ity of components of electromagnetic response measured 
with the receiver array. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention and its advantages may be more 
easily understood by reference to the following detailed 
description and the attached drawings in which: 

FIG. 1 is a perspective view of the preferred embodiment 
of the layout of the source and receiver apparatus used in the 
present invention for remote estimation of reservoir resis- 
tivities; 

FIG. 2 is a perspective view of an alternative embodiment 
of the layout of the source and receiver apparatus used in the 
present invention for remote estimation of reservoir resis- 
tivities; 

FIG. 3 is a plot of the axial vertical electric field complex 
magnitude |EJ as a function of the ratio of inner electrode 
radius over depth, a/d, for various values of the ratio of the 
inner electrode radius over electromagnetic skin depth, a/6; 

FIG. 4a is a plot of the axial vertical electric field complex 
magnitude |EJ versus the ratio of inner electrode radius over 
the depth, a/d; 
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FIG. 4b is a plot of the axial vertical electric field complex 
magnitude |Ej versus the ratio of inner electrode radius over 
the electromagnetic skin depth, a/6; 

FIG. 5a is a plot of the axial vertical electric field complex 
5 magnitude |Ej versus the ratio of outer electrode radius over 
the inner electrode radius, b/a; 

FIG. 5b is a plot of the total electrode current versus the 
ratio of outer electrode radius over the inner electrode 
radius, b/a; 

FIG. 6 is a perspective view showing the configuration of 
the alternative embodiment of the layout of the source and 
receiver apparatus used in the present invention for remote 
estimation of reservoir resistivities, as used in the example; 
15 FIG. 7 illustrates a bipolar square current waveform as 
output by the electromagnetic power source for use in the 
method of the present invention; 

FIG. 8 is a plot of the real part of the calculated radial 
component E r of the surface electric field response from the 
20 example; and 

FIG. 9 is a flowchart illustrating the processing steps of an 
embodiment of the method of the present invention for 
surface estimation of reservoir properties. 

While the invention will be described in connection with 
25 its preferred embodiments, it will be understood that the 
invention is not limited thereto. On the contrary, it is 
intended to cover all alternatives, modifications and equiva- 
lents that may be included within the spirit and scope of the 
3Q invention, as defined by the appended claims. 

DETAILED DESCRIPTION OF THE 
INVENTION 

The invention is a method whereby the average vertical 

35 and horizontal formation resistivities of a hydrocarbon res- 
ervoir are remotely mapped from the land surface or the 
seafloor, using low-frequency electromagnetic waves con- 
strained by seismic depth imaging and other a priori infor- 
mation about the subsurface. The invention overcomes the 
electromagnetic low vertical resolution problem by a com- 
bination of data acquisition and processing steps that are 
targeted at mapping the resistivity of a previously located or 
prospective reservoir. 

One embodiment of the layout of the apparatus used in the 

45 invention is shown in FIG. 1. In general, the invention uses 
the following four features that are synergistic in their 
combination: (1) a high-current multi-mode optimized elec- 
tromagnetic source, (2) a multi-component receiver array, 
(3) 3-D wave-equation data processing, and (4) reservoir 

50 properties estimation and mapping. These four features will 
be described in turn. 

In this embodiment of the present invention, two continu- 
ously grounded electrodes 4, 5, each consisting of one or 
more uninsulated electrical conductors, are buried at or 

55 within the near surface of the earth or the seafloor 1 in 
concentric circles of radii a and b respectively. Preferably, 
the electrodes are buried in the top zero to three meters. 
Alternatively, for offshore applications, the electrodes may 
be suspended or towed in the seawater above the seafloor as 

60 described by L. J. Srnka in U.S. Pat. No. 4,617,518. This 
electrode configuration provides for substantially maximiz- 
ing a transmitted vertical electrical current to a reservoir 
target 3 at depth d within the earth 2. In an alternative 
embodiment, an optional insulated circular loop 6, consist- 

65 ing of one or more electrically insulated conductors, is 
arranged at the earth's surface or seafloor 1 in a circle of 
radius c concentric to the two continuously grounded elec- 
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trodes 4, 5, for inducing horizontal electrical currents at the depth in meters. Here, p is average earth resistivity and f is 

reservoir 3. Although the insulated circular loop 6 is shown source excitation frequency. 

positioned between the two continuously grounded elec- FIGS. 4a and 46 are plots of the axial vertical electric field 

trodes 4, 5 in FIG. 1, the insulated circular loop 6 could be complex magnitude |EJ versus the ratio of inner electrode 

positioned inside the inner grounded electrode 4 or outside 5 fadius ovef the d h> a/d afld me fatk) of ^ electrode 

the outer grounded electrode 5. This insulated loop source 6 radius Qver ^ electromagnetic ^ depth) a/8> respectively, 

is used to augment natural background electnc and magneUc F , GS 3 and ^ show mat the axial | £j ^ maximized when 

field variations in the earth, to provide additional induced a/d _ 3/2 and ^ m Xhus> )E j ^ maximized when 

horizontal currents at the reservoir depth. The two grounded mM y^/A & shoW n in FIG. 46. 

electrodes 4, 5 and the optional insulated source loop 6are 10 ^ ¥ ^ 0 _ _ , , . . . . , i 

connected to one or more variable-frequency (preferably, ^ *° d 5 * show < he dependence of axial |Ej and 

lO^-lO 4 Hz) high current (preferably, lOMO^Amperes) total ^ ctIod * current on the ratio of outer electrode radius 

electrical power sources and controllers 7 by means of over the inner electrode radius b a, per unit source voltage, 

connecting cables 8, 9, preferably positioned radially. In the value of b/a = 9 maximizes |Ej at the reservoir, 

case of multiple sets of power sources and controllers 7, the 15 Preferably, the radius c of the insulated nng source 6 

connecting cables 8, 9 are preferably positioned equally (vertical magnetic dipole) is c^a, based on results from B. 

around the circumferences of the grounded electrodes 4, 5. R Spies (1989, op. cit.) for electromagnetic inductive 

The power sources and controllers 7 may be located at the sounding in the near-field zone of a vertical magnetic dipole 

land surface or seafloor 1. Alternatively, in the case of source. 

offshore surveys, the power sources and controllers 7 may be 2 o An alternative embodiment of the grounded electrodes is 
located at the sea surface, or within the body of the sea. The shown in FIG. 2. The two concentric ring electrodes 4, 5 of 
power sources and controllers 7 provide for selective exci- respective radii a and b are replaced by six or more linear 
tation of the grounded electrodes 4, 5 and insulated wire grounded electrodes 11 of equal lengths L«b-a. Preferably, 
loop 6, modifying the frequency as required to maximize the these electrodes 11 will be continuously grounded to the 
subsurface response. 25 eartn along their entire individual lengths L. Alternatively, 
The dimensions and probing frequency for a given res- the electrodes 11 may be only partially grounded, that is, 
ervoir depth and average earth resistivity, plus the corre- continuously grounded only within some distance y=L/2 as 
sponding electrical impedance of the grounded electrodes 4, measured from the radii a and b, respectively, as shown in 
5, are calculated by numerically solving the uninsulated FIG. 6. Preferably, the electrodes 11 are placed along radii 
buried low-frequency electromagnetic antenna problem. 30 separated by equal angles of not more than 60°, whose 
Preferably, this problem is solved using the methods of R. inward radial projections intersect at the center of the 
W. P. King and G. S. Smith, Antennas in Matter, MIT Press, electrode array. Preferably, each linear electrode 11 is con- 
Cambridge, 1981. Preferably, the solution is implemented nected at each of its ends (r-a and r=b) to a continuously 
by means of a 3-D frequency-domain computer program in grounded linear terminating electrode 12 that is substantially 
which the surface potentials, current densities, and electric 35 orthogonal (preferably, 90°± 10°) to the connected radial 
fields are found by solving Maxwell's equations using 2-D electrode 11. Preferably, the length of each terminating 
complex Fourier transforms at each depth interval. Bound- electrode 12 is not more than L/10. Preferably, one or more 
ary conditions for solving the problem are applied to enforce power sources and controllers 7 are connected to the radial 
the condition that, at each frequency, the total current linear electrodes 11 near the midpoints L/2 of the electrodes 
leaving the inner ring 4 equals the total current captured by 40 11 within a distance of ±1710. If more than one power source 
the outer ring 5, and that the voltage difference between the 7 is used simultaneously, the multiple sources 7 operate in 
rings 4, 5 is conserved. Preferably, the values of radii a and a synchronized manner to supply electrical current to each 
b are determined by substantially maximizing the vertical electrode 11. Preferably, source synchronization is such that 
and horizontal electric fields at the symmetry axis of the the total phase variations of the six or more source currents 
concentric rings 4, 5 (radius r=0) at the depth of the center 45 do not exceed 0.1 degree and the total amplitude variations 
of the reservoir. FIGS. 3, 4a, 4b, 5a, and 56 show results of the source currents do not exceed 0.1 percent. In this 
from this calculation for sample input values. The sample alternative embodiment employing grounded electrodes 11, 
input values affect the absolute values of the electric fields the optional insulated circular wire loop 6 may also be used, 
and currents, but do not affect the dimensionless scaling as described above and shown in FIG. 6. The power sources 
parameters used to optimize the source electrodes. Thus, in 50 ? operate in a discrete-frequency ("frequency domain") or a 
FIGS. 3, 4a and 46, the shape of the curves would remain the variable-sequence alternating wave ("time domain") tran- 
same for the same ratio of outer electrode radius to inner sient manner. In both cases, the polarity of the source 
electrode radius b/a-8, while the absolute magnitude of the currents is reversed periodically (preferably, 10" 4 to 10 4 
curves would change for different values of average earth seconds) as in standard commercial practice well known to 
resistivity and source excitation voltage. Similarly, in FIGS. 55 tDOSe of ski11 m the **U m order t0 minimize electrode 
5a and 56, the shape of the curves would remain the same polarization effects. 

for the same ratio of depth to inner electrode ratio d/a=2/3, Electromagnetic responses are collected by an array of 

while the absolute magnitude of the curves would change for multi-component receivers 10 positioned at the surface of 

different values of average earth resistivity and source the earth or at the seafloor 1, as shown in FIGS. 1 and 6. 

excitation voltage. These ratios and values were selected for 60 Preferably, two orthogonal horizontal electric fields, two 

illustrative purposes only and should not be taken as limi- orthogonal horizontal magnetic fields, and a vertical mag- 

tations of the method of the present invention. netic field are measured, when the array of receivers 10 is 

FIG. 3 gives the axial (radius r«0) vertical electric field positioned on land. Preferably, the orthogonal horizontal 

complex magnitude |Ej (in millivolts/meter) as a function of directions are the same for all receivers 10, Preferably, an 

a/d for various values of the ratio a/6, per unit voltage (1 65 additional vertical electric field component is also measured 

Volt) of source excitation between the grounded electrode when the array of receivers 10 is positioned at the seafloor 

rings 4, 5, where 5- 503x(p/f) 1/2 is the electromagnetic skin 1. Preferably, receiver signal amplitudes and phases are 
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measured with an accuracy greater than or equal to 0.1 %, and from many source repetitions, combined with local 
relative to the source fields, using commercially available noise measurements and signal cross-correlation techniques, 
broadband (preferably, 10" 4 to 10 4 Hz) electric and magnetic are used within these standard methods to achieve noise 
sensors that have standard high sensitivity, and a receiver suppression. Preferably, such techniques are applied to the 
system dynamic range given by a capacity of 24 bits or 5 data to produce a signal-to-noise ratio greater than or equal 
more. Preferably, the magnetic field sensors have phase to 1 and signal accuracy greater than or equal to 1.0% for 
accuracy greater than or equal to 0.1 degree over the each electromagnetic component used within the multi- 
frequency range used for the survey. In the preferred component data inversion. 

embodiment of the method of this invention, these multi- The electromagnetic source signature (source-generated 

component responses at each receiver 10 are also recorded JQ ^ suppressed automatically by the self-canceling field 

when the grounded and ungrounded sources 7 are turned off, geometry of the grounded electrodes (4 and 5 in FIG. 1, 11 

to measure the earth's electromagnetic response to the an d 12 in FIG. 2), as described by Mogilatov and Balashov 

natural background electric and magnetic fluctuations and (1996, 0 p. cit.). Additional suppression (deconvolution) of 

also to measure the electromagnetic noise environment. source effects is accomplished by normalizing the receiver 

Preferably, electromagnetic responses are measured at 15 data to the background earth response using standard indus- 

each receiver site over a grid having receiver spacing try techniques well known to one of skill in the art, such as 

intervals x and y^0.5 d, where d is the vertical distance described in M. Zhdanov and G. Keller, (1994, op. cit.) or 

(depth) from the land surface or seafloor 1 to the reservoir K. M. Strack (1992, op. cit.). Alternatively, this suppression 

3, as shown in FIGS. 1, 2, and 6. The x and y intervals may may be accomplished by normalizing (cross-referencing) 

differ. Alternatively, linear receiver arrays (one or more 20 the data using apparent resistivity functions for a layered 

parallel lines of receivers 10) may also be used. The linear earth as described in T. G. Caldwell and H. M. Bibby, 

arrays may also be arranged in swath geometry, in which Geophys. 7. Int., v. 135, 817-834, 1998. 

case the receiver data may be summed in the cross-line Preferably, the data from the array of receivers 10 are 

direction. Preferably, receivers 10 are positioned in a grid converted ("inverted") from time or frequency domain elec- 

pattern, as described above, over the entire area from the ^ tromagnetic responses into a 3-D resistivity depth image of 

center of the source array out to a radial distance r=b. The the earth by the application of iterative 3-D model-guided 

few receivers 10 shown in FIGS. 1 and 6 illustrate the nonlinear electromagnetic inversion methods that incorpo- 

inter-receiver spacing dimensions without showing the rate geometrical and electrical parameter constraints, as will 

extent of the preferred coverage, and thus should not be be described below. Finite-difference and finite-element 3D 

taken as a limitation of the present invention. This position- 30 models may be used. Inversion methods used in this inven- 

ing maximizes the data sensitivity to the reservoir 3 and tion include standard techniques such as quasi-linear regu- 

other electrical structures within the earth 2 near the reser- larized methods, such as described in M. S. Zhdanov and S. 

voir 3, and provides the greatest depth of penetration (the Fang, Radio Sci., v. 31, 741-754, 1993 and fully nonlinear 

"near-field" response) for a given electrical or magnetic conjugate gradient or Gauss-Newton methods, such as 

dipole source moment M. 35 described in G. A. Newman and D. L. Alumbaugh, Report 

The spatial positions and orientations at the earth's sur- SAND96-0582, Sandia National Laboratories, 1996; and 
face or seafloor 1, of the source electrodes 11, and of the Geophys. J. Int., v. 128, 345-354, 1997; Alumbaugh and 
receivers 10 in the array are measured. Preferably, the Newman op. cit., 1997 and SEG Annual Meeting Expanded 
positioning and orientation is accomplished using geodetic Abstracts, v. 1, 456-459, 1998. Preferably, the positions and 
methods standard in the industry and well known to those of 40 strengths of all source currents applied at or within the 
skill in the art. These geodetic methods may include differ- surface of the land surface and at, above, or within the 
ential and kinematic GPS (Global Positioning Satellite), and seafloor 1 are explicitly included in the inversion, by means 
acoustic transponders in an offshore application. Preferably, of Green's functions or other standard mathematical tech- 
maximum allowed position uncertainties are ±0.001 d in the niques that are well known to one of skill in the art. 
vertical and the two horizontal directions. Preferably, maxi- 45 Preferably, separate inversions are performed for receiver 
mum allowed orientation uncertainties are ±0.10 degrees in data collected using the grounded electrode sources (4 and 
the vertical and in the two horizontal orientations. 5 in FIG. 1, 11 and 12 in FIG. 2) and for receiver data 

Preferably, the multi -component electromagnetic receiver collected using the insulated loop source 6 and the earth's 

data are processed using full wave-equation methods. This natural background magnetotelluric source when the other 

3-D processing includes, but is not limited to, data noise 50 sources are turned off. Preferably, joint inversions of 

suppression, source deconvolution, and model-guided inver- receiver data collected using any combination of the 

sion. Both frequency-domain and time-domain methods are grounded source (4 and 5 in FIG. 1, 11 and 12 in FIG. 2), 

used, depending upon the methods used for the data acqui- insulated source 6, and magnetotelluric source are also 

sition. Alternatively, electromagnetic wave-equation migra- performed. 

tion may be used such as that described by M. Zhdanov et 55 Preferably, the spatial positions of geometrical constraints 
al, Exploration Geophysics, v 26, 186-194, 1995; M. are obtained from surfaces, such as horizons and faults, 
Zhdanov and O. Portniaguine, Geophys. J. Int., v 131, interpreted in dense 2-D or in 3-D depth-converted stacked 
293-309, 1997; and M. Zhdanov et al, SEG Annual Meeting seismic reflection data. Preferably, standard industry seismic 
Expanded Abstracts, v. 1,461^168,1998. interpretation packages, such as Geoquest IESX©, Para- 
Preferably, standard electromagnetic industry data pro- 60 digm GeoDepth©, or Jason Workbench©, are used to pro- 
cessing techniques such as those described by M. N. Nabig- duce the interpreted seismic surfaces, to tie the seismic depth 
hian (1988, op. cit.); K.-M. Strack (1992, op. cit.); G. Buselli data to well log, gravity, magnetic, and other geoscience 
and M. Cameron, Geophysics, v. 61, 1633-1646, 1996; and data, and to transfer these depth surfaces to the 3-D elec- 
G. D. Egbert Geophys. J. Int., v. 130, 475^496, 1997 are tromagnetic inversion starting model. Resistivity values for 
used forsuppressionofboth natural background and human- 65 the initial electromagnetic depth model for geologic units 
generated electromagnetic noise. Preferably, data redun- bounded by the interpreted seismic surfaces are produced by 
dancy from multiple-receiver multi-component responses any of a number of standard industry methods well known 



US 6,603,313 Bl 

11 12 

to one skilled in the art. These methods include ties to log hydrocarbon pore volume, net-to-gross, water saturation, 
data; extrapolation from regional data bases; application of and other reservoir properties from the and 
empirical resistivity transforms using seismic intervals, well inversion measurements. The derivation uses the facies 
sonic velocities, or acoustic impedances; and initial layered- model of rock properties distributions combined with Arch- 
earth (1-D) resistivity inversion derived from the collected 5 ie's Equations for the electrical resistivity of a porous rock 
electromagnetic receiver data. Preferably, constraints are containing fluid in the pore spaces, relative to values 
enforced during the inversion using standard industry within the same geologic unit outside of the reservoir, 
techniques, such as described in M. A. Meju, Geophysical The invention described above is designed to provide an 
Data Analysis: Understanding Inverse Problems and order of magnitude improvement in subsurface vertical 
Theory, Society of Exploration Geophysicists, 1994. These 10 electromagnetic resolution over current technology, 
standard techniques include Tikhanov regularization, Baye- FIG. 9 is a flowchart that illustrates a preferred embodi- 
sian methods, sharp-boundary approaches (G. Hoversten et ment of the method of the invention for surface estimation 
al, 1998, op. cit.), equivalent integral conductance and of reservoir properties of a subsurface geologic formation, as 
resistance methods, and minimum gradient support tech- just described. First, at step 900, location of the subsurface 
niques (O. Portniaguine and M. Zhdanov, 1998, op. cit.). 15 geologic formation is determined, using geological and 

In an alternative embodiment of the invention, interpre- geophysical data in the vicinity of the subsurface geologic 
tation of the inverted resistivity depth cubes ("inversions") formation. Next, at step 902, average earth resistivities 
includes comparison of the 3-D resistivity-depth values with above, below, and horizontally adjacent to the subsurface 
interpreted 3-D seismic features and all mapped attributes geologic formation is determined, using geological and 
derived from the seismic data (pre-stack and post-stack). 20 geophysical data in the vicinity of the subsurface geologic 
Preferably, the separate and mathematically joint electro- formation. Next, at step 904, dimensions for a high-current 
magnetic inversions produced from the grounded electrode multi-mode electromagnetic source are determined to sub- 
sources and from the insulated loop source and magnetotel- stantially maximize transmitted vertical and horizontal elec- 
luric source are compared and combined with each other and trie currents at the subsurface geologic formation, using the 
with the seismic reflection features and seismic attributes to 25 location and the average earth resistivities. Preferably, the 
estimate the fluid type, hydrocarbon pore volume, dimensions are calculated by numerically solving the unin- 
saturation, and the shahness factor (net pay-to-gross reser- sulated buried low-frequency electromagnetic antenna 
voir thickness ratio) within the reservoir over its mapped problem, as described previously. Next, at step 906, probing 
extent. The reservoir may be seismically defined by a frequency for a high-current multi-mode electromagnetic 
combination of stratigraphic or structural closure or limits of 30 source is determined to substantially maximize transmitted 
mapped seismic attributes. vertical and horizontal electric currents at the subsurface 

The preferred method of this alternative embodiment of geologic formation, using the location and the average earth 
the invention to estimate fluid type, hydrocarbon pore resistivity Again, the probing frequency preferably is 
volume, saturation and shaliness factor is as follows. The calculated by numerically solving the uninsulated buned 
resistivity inversion at the reservoir depth produced from 35 low-frequency electromagnetic antenna problem as 
electromagnetic receiver data collected using the grounded de f cr ' bed previously Alternatively, iterated 3-D modeling 
electrode sources is used to measure the vertically averaged calculations of the subsurface geologic formation s electro- 
resistivity p v „, within the reservoir 3. The resistivity inver- magnetic response may be used to verify the dimensions and 
sion at the reservoir depth produced from electromagnetic P robm g frequency of the high-current mulu-mode electro- 
receiver data collected using the insulated source or the 40 magneUc source m steps 904 and 906. Next, at step 908 the 
magnetotelluric source is used to measure the horizontally electromagnetic source is activated at or near the surface, 
averaged resistivity p horiz within the reservoir. A facies approximately centered above the subsurface geologic for- 
model of the reservoir is derived from the seismic matl0n - Next '. at ste P 910 > a Polity of components of 
interpretation, geologic concepts, and available facies data electromagnetic response are measured with a receiver 
(such as well logs and data bases). This facies model is 45 array Preferably, when the array of receivers 10 ,s posi- 
combined with the PverI and Pjtarii[ inversion measurements tto f d 00 a ° d > *"° orthogonal horizontal electric fields, two 
to estimate products and ratios of the permeable bed resis- orthogonal horizontal magnetic fields, and a vertical mag- 
tivity p„, the summed permeable bed thickness divided by netic field are measured Alternately, when the array of 
the total reservoir interval thickness ntg ("net-to-gross"), receivers 10 is positioned offshore, an additional vertical 
and the impermeable bed resistivity p sll . For a reservoir 50 e ectnc field is measured. Next, at step 912, geometneal and 
facies model comprised of uniform values p„ for the per- electrical parameter constraints are determined, using the 
meable beds, and a different but uniform value Pj „ for the geological and geophysical data. Next, at step 914, the 
impermeable beds, then, as it is known in the industry: electromagnetic response is processed using the geometrical 

and electrical parameter constraints to produce inverted 

P^,j-P,^"'«+P,*x(i-"'8) (t) 55 vertical and horizontal resistivity depth images. Preferably, 

, . m me components of the electromagnetic response are pro- 

cessed using full 3-D wave-equation methods, as described 

Equations (1) and (2) contain three unknown averaged previously. 1-D inversion of the electromagnetic response is 

reservoir parameters: p„, p sh , and ntg. Estimates of p sh used to verify the average earth resistivities above, below, 

within the reservoir interval, derived independently from the 60 and horizontally adjacent to the subsurface geologic 

facies model or available facies data, are used next to derive formation, as determined in step 902. Finally at step 916, the 

and map the two remaining unknown values p„ and ntg over inverted resistivity depth images are combined with the 

the spatial extent of the reservoir. Reservoir fluid type, geological and geophysical data to estimate the reservoir 

hydrocarbon pore volume or water saturation are then properties. Details of the preferred method of inversion are 

derived from the mapped p„ value within the area of the 65 described later in conjunction with the following example, 

seismically defined reservoir. Statistical methods including The following example illustrates the application of the 

Monte Carlo inversions may also be used for deriving invention for onshore (land) hydrocarbon reservoir resistiv- 



US 6,603,313 Bl 

13 14 

ity mapping. After 3-D seismic data in the survey area are Nine power sources/controllers 7 are placed at the surface 
acquired, interpreted, and converted to the depth domain, the of the earth 1. Each source/controller is powered by con- 
prospective reservoir is identified (depth d and extent 1). nection to a municipal power grid, if available, or is powered 
Knowledge of the earth's electrical resistivity for the survey by one or more generators in the field survey area. Each 
area, averaged over intervals of O.lOxd, from the earth's 5 source/controller is nominally rated at 300 kVA, with out- 
surface to three times the reservoir depth (3xd) and five P uts of 120 VAC and 2500 A (rms). One source/controller is 
times the reservoir extent (5x1), is gathered using existing situated at any position along the circumference of the 
electromagnetic survey data and well logs, or is estimated ™^* x ? loo P 6 ' an f d * connected by a coaxial power 
usinggeologicbasinanalogs.Thediametersofthegrounded cable 9 at the surfaoe of u the earth 10 ^ insulated loop 
electrodes are calculated by numerically solving The unin- 10 The pmaimng eight power sources/controllers 7 are 

. A , . . - . - , ^ ' • f . placed within a distance of L/10 of the midpoints of the 
sulatedbuned low-frequency electroma^euc antenna prob- * ^ M electrodes, as shown in FIG. 2. These 
lem i as discussed above, or by -iterated 3D electromagnetic ^ f wtctstconKoiim are connected to the 
modeling, using the reservoir depth and vertically averaged groimded electrodes by means of coaxial or single- 
layered-earth resistivities as inputs. The diameter of the cond uctor power cables 8. The satellite Global Positioning 
optional insulated loop electrode is determined using stan- 15 System (GPS) signal is used to monitor and synchronize the 
dard methods known in the art. phases of all the sources. Alternatively, if the number of 

FIG. 6 shows land source and receiver configurations for p0 wer sources/controllers 7 is limited, or if survey logistics 

a target reservoir 3 identified seismically at d=1000 meters or terrain difficulties make simultaneous use of the eight 

depth to top of reservoir, having an average lateral extent radial grounded electrode positions impractical or too costly, 

(radius) 1/2=1250 meters. Eight (8) partially grounded radial 20 one or more power sources/controllers may be used to 

electrodes U and connected terminating electrodes 12, as energize separately the eight partially grounded and the 

described also in FIG. 2, are deployed in a radial array in insulated loop source, in any sequential order, 

conjunction with an insulated loop source 6. The geometri- The partially grounded electrodes 11 and the terminating 

cal center of the grounded electrode array (intersection of electrodes 12 each consist of three uninsulated size 4/0 

their 8 radius lines) and the center of the insulated loop are 25 multi-strand copper wires. The grounded wires that com- 

positioned at the surface of the earth 1 vertically above the prise the grounded portion of each partially grounded elec- 

center of the reservoir target. The grounded electrodes are trode element 11 are buried in parallel within the top 1.0 

positioned symmetrically around the circumference of the meter of the earth's surface by means of manual digging or 

source array, each separated by an angle of 45±1 degrees standard mechanical cable-laying devices. The ungrounded 

from the adjacent electrode as measured from the center of 30 portions of each of the radial partially grounded electrodes 

the source array. The source dimensions are a«1500 meters, 11 consist of three uninsulated size 4/0 multi-strand copper 

b=6000 meters, y=90 meters, and c=1000 meters. The value wires that are connected to the uninsulated buried electrode 

Y is determined from the calculation of vertical current wires comprising the grounded portions. The ungrounded 

leakage from a continuously grounded bipole antenna of portions of each radial electrode are laid on the surface of the 

length L, using the method described above to numerically 35 earth. Electrical contact of the grounded radial electrodes 

solve the uninsulated buried low-frequency electromagnetic and the terminating electrodes is maintained with the earth 

antenna problem. This shows that most of the current leaves by periodically wetting the buried electrode areas with 

the grounded wire within a distance iL/5 at each end of the water, as needed according to local ground moisture condi- 

antenna. The grounded terminating electrodes 12 each have tions. The loop source 6 consists of one single-conductor 

a length of 30 meters. The grounded array and the insulated 40 multi-strand insulated size 4/0 copper wire. Power connec- 

loop are not moved during the survey. Alternatively, if the tion cables 8 and 9 are electrically rated according to U.S. 

number of power sources/controllers 7 is limited, or if NEMA (National Electrical Manufacturing Association) 

survey logistics or terrain difficulties make simultaneous use codes and standards to carry the current delivered to the 

of the eight radial grounded electrode positions impractical grounded electrodes 11,12 and to the insulated loop 6, 

or too costly, the eight radial partially grounded electrode 45 respectively. 

positions are occupied sequentially in groups of one or more Electromagnetic receivers 10, such as Electromagnetic 

positions, in any sequential order. Instruments, Inc. (EMI) type MT-24/NS™ or equivalent, are 

A preferred procedure is to obtain substantially optimal positioned over the surface of the earth 1 within a radial 

parameter values to substantially maximize the electric field distance r=(x 2 +y 2 ) 1/2 =5000 meters from the center of the 

at the reservoir depth. However, as an alternative procedure, 50 array, but not within 25 meters of any grounded electrode 11, 

a sub-optimal aspect ratio b/a could be used to reduce 12 or the insulated loop 6, to minimize source-generated 

electrode cost, installation effort, and survey permitting. For noise and saturation of the receiver signals. The receivers are 

instance, an aspect ratio b/a =4 could be used. Use of this positioned on a uniform grid as shown in FIG. 6, with a 

value for b/a would result in a 24.5% reduction in vertical lateral spacing of x«y=100 meters, within a radius of 2000 

electric field at the reservoir target, as shown in FIG. 5a, and 55 meters from the center of the array, and on a uniform grid 

a corresponding reduction in the electromagnetic responses with a lateral spacing of x-y=300 meters from a radius of 

of the reservoir to the grounded electrode excitation as 2500 meters to a radius of 5000 meters from the center of the 

measured at the surface receiver array 10. array. Each five-channel receiver measures two components 

Assume a vertically averaged resistivity of the earth of (x and y directions) of the horizontal electric field, two 

value p^-1 Ohm-m. Then the central operating frequency of 60 components (x and y directions) of the horizontal magnetic 

the grounded electrode array is derived from d/5=9/4 and field, and one component (z direction) of the vertical mag- 

d=2250 meters, or f=0.050 Hz. The output bandwidth of the netic field. The receivers are modified by standard industry 

grounded electrode sources is 0.005 ^f^5.0 Hz. Using the methods including feedback stabilization so that the phase 

analysis of B. R. Spies (1989, op. cit.), the central operating accuracy of the magnetic field induction sensors (EMI type 

frequency of the insulated loop source is set by d/5-1, or 65 BF-4™ or equivalent) is greater than or equal to 0.10 

f=0.253 Hz. The output bandwidth of the insulated loop degrees in the full frequency range of the survey 

source is 0.025 if i 25 Hz. (0.005 if i25 Hz). 
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The five-component receivers are deployed simulta- frequency-wavenumber domain using standard industry 2-D 
neously in large groups (16 or more) within the survey area, Fourier and Radon transform techniques. The first set of 
with as many receiver groups deployed as possible and magnetotelluric data and the second set of vertical magnetic 
practical for the local conditions of the survey (e.g. terrain dipole data are merged together in the frequency- 
difEculties, logistical support). Data are gathered for each 5 wavenumber domain, for each electromagnetic tensor corn- 
receiver group by a central processing unit (EMI type ponent of the data. The merged magnetotelluric and vertical 
FAM/CSU™ or equivalent). Differential GPS geodetic magnetic dipole data sets are inverted, and the grounded 
methods are used to measure the positions (x, y, z) of all radial electrode data set is inverted separately. Then the 
receivers to within 0.1 meters accuracy. The GPS signal is merged magnetotelluric and vertical magnetic dipole data 
also used for phase synchronization (timing) of all receiver 10 and the grounded radial electrode data are inverted jointly, 
data. as discussed in D. Jupp and K. Vozoff, Geophys. 

The receiver data are collected in three ways. First, the Prospecting^ v. 25, 460-470, 1977. The magnetotelluric 
receiver data are collected as time records with all sources data, the vertical magnetic dipole data, and the grounded 
7 turned off, to record zero excitation currents. These data radial electrode data are also inverted separately. All data 
are collected over a length of time that is sufficient to record 15 inversions use the 3-D frequency-domain finite-difference 
raw stacked magnetotelluric data having three-sigma errors fully nonlinear methods of G. A. Newman and D. L Alum- 
less than or equal to 5% over the frequency range baugh (1996, 1997, op. cit.), modified to allow for the 
0.0025 ^f ^25 Hz. Typically, collection of this data will take geometries of the grounded radial electrode and the insu- 
1-10 days, depending upon local conditions and the logistics lated loop source current arrays. Depth and parameter value 
of receiver deployment. This first set of receiver data is 20 constraints are enforced during the inversion, using sharp- 
magnetotelluric data. Second, the insulated loop source is boundary methods (G. Hoversten et al, 1998, op. cit.) and 
energized using a standard electromagnetic industry bipolar integral resistance and conductance bounds within the 
square wave current from its attached power source/ update region of the nonlinear inversion 3-D mesh that 
controller 7, as shown in FIG. 7. In this preferred embodi- contains the reservoir target, combined with minimum- 
ment of the method of the present invention, current pulse 25 gradient support techniques (O. Portniaguine and M. 
on -time, Tj equals current waveform off-time, T 2 , that is, Zhdanov, 1998, op. cit.). The nonlinear inversion update 
Tj-T 2 , but this is not a limitation of the method. Other region is centered on the target reservoir, and extends 100 
source current waveforms may also be used for the insulated meters above and below the reservoir and 200 meters 
loop source current, including sinusoidal waveform combi- laterally from each reservoir edge. 

nations and pseudo-random sequences as well known to one 30 The starting model for both the merged magnetotelluric 
skilled in the art, provided the insulated loop source fre- and vertical magnetic dipole data inversion and the 
quency range is as stated. The receiver responses are col- grounded radial electrode data inversion is an interpreted 
lected using time-domain measurements acquired during the seismic depth model in which the mechanical properties 
current waveform off-time, T 2 in FIG. 7. The duration of the (primarily the interval acoustic impedances) are replaced 
current pulse on-time, T a in FIG. 7 (and hence also the 35 with resistivity estimates. The resistivity estimates may 
off-time T 2 ), is set at three values, 0.01, 1.0, and 10.0 come from electromagnetic survey data, well logs, empirical 
seconds. Sufficient repetitions (typically 50 to 1000) of the relations to seismic parameters, or geologic basin analogs, as 
loop source current are made at each on-time value so that described above. The inversions are performed by means of 
the raw stacked data time series data have three-sigma errors a digital electronic computer of the massively parallel pro- 
less than or equal to 1% over the frequency range 40 cessor (MPP) type, or alternatively using a network of 
0.025 ^f ^25 Hz. This second set of receiver data is vertical electronic digital computers that mimic an MPP computer, 
magnetic dipole data. Third, the insulated loop source is After the separate magnetotelluric, vertical magnetic dipole, 
turned off (zero current) and the eight partially grounded and grounded radial electrode data inversions are completed, 
electrodes 11, 12 are simultaneously energized in phase. the magnetotelluric — vertical magnetic dipole and grounded 
Alternatively, if the number of power sources/controllers is 45 radial electrode data are inverted jointly. The five respective 
limited, the partially grounded electrodes are energized 3-D depth cubes of inverted resistivity (magnetotelluric, 
separately in groups of one or more, in any sequential order. vertical magnetic dipole, grounded radial electrode, merged 
Each of the eight power sources/controllers 7 produces a magnetotelluric — vertical magnetic dipole, and merged 
standard electromagnetic industry bipolar square wave cur- magnetotelluric — vertical magnetic dipole — grounded 
rent pulse, as shown in FIG. 7, with the duration of the 50 radial electrode) are compared, and the ratios of their 
current pulse on-time, T 1 in FIG. 7, (and hence T 2 ) set at resistivity values are formed at each depth location using 3D 
three values, 0.05, 5.0, and 50.0 seconds. Other source visualization methods. Finally, values of p„, p shy and ntg are 
current waveforms may also be used for the grounded source derived for the reservoir interval using the methods 
current, including sinusoidal waveform combinations and described above, and are mapped. These mapped values are 
pseudo-random sequences as well known to one skilled in 55 interpreted in conjunction with the 3-D seismic data and its 
the art, provided the grounded electrode source frequency attributes. 

range is as stated. Sufficient repetitions (typically 50 to FIG. 8 shows the complex magnitude of the calculated 
1000) of the grounded electrodes' source currents are made radial component E r =(E r 2 +E y 2 ) 1/2 of the surface electric 
at each on -time value so that the raw stacked data time series field response from the example target reservoir described 
data have three-sigma errors less than or equal to 1% over 60 above, due to excitation by the grounded electrode array, 
the frequency range 0.005 ^f ^5 Hz. This third set of The example reservoir is assumed to have a vertical thick- 
receiver data is grounded radial electrode data. ness of 20 meters and a vertically averaged resistivity of 100 
The three sets of receiver data are processed in the Ohm-m. The electromagnetic response was calculated using 
following way. After noise suppression using standard the SYS EM 3-D integral equation computer code developed 
industry methods as described above, the second set of 65 at the University of Utah's Consortium for Electromagnetic 
vertical magnetic dipole data and the third set of grounded Modeling and Inversion. This electric field component 
radial electrode measurements are converted to the complex response is normalized to the uniform earth (halfspace) 



